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ABSTRACT 

Measurements of the asymmetry of the emission peaks in the core of the Ca II 
H line for 105 giant stars are reported. The asymmetry is quantified with the 
parameter V j R, defined as the ratio between the maximum number of counts in 
the blueward peak and the redward peak of the emission profile. The Ca II H and 
K emission lines probe the differential motion of certain chromospheric layers in 
the stellar atmosphere. Data on V j R for the Ca II K line are drawn from previous 
papers and compared to the analogous H line ratio, the H and K spectra being 
from the same sets of observations. It is found that the H line (V/R)h value is 
+0.04 larger, on average, than the equivalent K line ratio, however, the difference 
varies with B — V color. Red giants cooler than B — V = 1.2 are more likely 
to have (V/R)h > (V/R)k, whereas the opposite is true for giants hotter than 
B — V = 1.2. The differences between the Ca II H and K line asymmetries could 
be caused by the layers of chromospheric material from which these emission 
features arise moving with different velocities in an expanding outflow. 



1. Introduction 

The chromospheric emission component of the Ca II K line has been studied extensively 
in the optical spectra of giant stars (e.g., Wilson 1976, Blanco et al. 1976, Stencel 1978, 
Kelch et al. 1978, Gray 1980, Middelkoop & Zwaan 1981, Middelkoop 1982, Rutten 1984). 
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The Ca II H line has not been used as often to probe the stellar atmosphere, though both 
lines exhibit emission in their cores that typically has a double-peaked profile. In the H 
line, the violet- and red-side emission peaks can be labeled H 2 y and H 2 r respectively. A re- 
absorption feature denoted H 3 is located at the core of the emission profile, and is produced 
higher in the chromosphere. The H 2 y and H 2 # peaks often do not have equal strength. To 
quantify this asymmetry, it is convenient to use the V/R parameter of Wilson (1976) and 
Stencel (1978), which we define as the ratio of intensities in the H 2 y and H 2 ^ emission peaks. 
Asymmetric emission profiles with V / R values less than unity can be interpreted as evidence 
for an outflow having a differential velocity field in the chromospheric region, with the H 3 
absorption feature being blueshifted with respect to the emission profile (Stencel 1978). It 
is possible that the Ca II H line probes a different level of the stellar chromosphere than 
the K line, depending upon the mechanism for setting the relative populations in the upper 
states of these transitions (e.g., Avrett 1966; Linsky 1970). With the greater opacity of the 
K line leading to formation at higher altitudes, these emission lines can potentially be used 
to probe the onset of mass outflows within two different layers of the chromosphere. 

In this paper we report the results from a spectroscopic program to measure the V/R 
parameter of the Ca II H line for a sample of red giant stars. The asymmetry of the Ca II 
H emission profile is compared to that of the K line. Our study uses spectra of 105 red 
giants observed during a ten-year period, and includes 98 stars previously investigated by 
Smith & Shetrone (2000, 2004) and Shetrone et al. (2008). The sample covers a range of 
stellar colors and luminosities in an effort to map out the behavior of the H emission as a 
function of location in the H-R diagram. In addition, multiple observations were obtained 
for a number of stars. Variability of the K emission line profile based on our data set has 
been discussed by Shetrone et al. (2008). Although time variability is not the intended 
subject of this paper, the individual observations listed in Table 1 allow an interested reader 
to judge the degree to which the H line asymmetry can vary on a timescale of months. 



2. Observations and Reduction 

High resolution spectra of the Ca II H and K lines for a sample of 105 field giants were 
obtained at McDonald Observatory. The spectra were acquired during several observing 
runs between October 1998 and June 2008 using either the Sandiford Cassegrain Echelle 
(CE) spectrometer on the 2.1-m Otto Struve Telescope or the 2d-coude spectrometer on the 
2.7-m Harlan J. Smith Telescope. For the 2.7-m we used a 1.2" slit which yielded a 2.0 
pixel resolution of R = 60,000. For the 2.1-m we used a 1.1" slit which yielded a 2.0 pixel 
resolution of R = 60, 000. This paper examines the Ca II H line found in a different order 
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of the same echelle spectra for which K line observations have been reported upon by Smith 
& Shetrone (2000, 2004) and Shetrone et al. (2008). To this material we add more recent H 
and K line data from the spectra obtained in June 2008. 

The observing procedure and data reduction process for the spectra obtained during 
observing runs from October 1998 through August 2007 are described in Smith & Shetrone 
(2000, 2004) and Shetrone et al. (2008). The most recent data were taken during the nights 
of 2008 June 22 and 23 using the 2d-coude spectrometer on the 2.7-m telescope. For this 
later run, the observing procedure was to take at least three exposures of each object with 
the maximum exposure time being 1200 s. Bias frames and flat field frames were taken 
at the end of each night. In the case of the 2.1-m program, since the spectrometer is a 
cassegrain instrument, exposures of ThAr arcs were made before and after observation of 
each star. By contrast, the 2d-coude spectrometer of the 2.7-m telescope is very stable, 
and for the observing runs with this instrument, ThAr arcs were generally obtained only 
at the start and end of each night. The spectra were bias corrected, flat-field corrected, 
and wavelength calibrated using routines within the IRAF echelle package. The individual 
spectra for each object taken on a given night were combined to remove cosmic rays. During 
the data reduction process, the number of observed electrons was preserved to keep the 
Poisson statistics valid for all coadded spectra. 

Using the IRAF fxcor and dopcor packages, the spectra were shifted to the velocity 
frame of the Hinkle spectral Atlas of Arcturus (Hinkle et al. 2000). The spectra were then 
divided by the Arcturus Atlas to obtain a continuum residual, which was fit by a low order 
Legendre polynomial, excluding the core of the Ca II H line in the fitted sample. The 
combined spectra were divided by the continuum fit to obtain normalized spectra. This 
normalization procedure matches the procedures found in Smith & Shetrone (2004) and 
Shetrone et al. (2008). The original spectrum, before continuum normalization, was also 
saved and referenced to determine the Poisson noise. 



3. Results and Discussion 

3.1. The Asymmetry of the Call H 2 Emission 

Our sample is comprised of 105 giants that show double-peaked Ca II H and K emission 
lines in their spectra. The stars observed are listed in Table 1 along with each date on which 
they were observed and the (V/R)h value for the H emission feature. The error (denoted 
en) on each (V/R)h value is calculated from the relation e 2 H = cry/Ij I2R + Ih2V (T r/ ^H2Ri as 
described in Shetrone et al. (2008). In this equation, the recorded photon counts in the H 2 y 
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and B.2R peaks of the unnormalized spectra are represented by Imv and Ih2R, respectively. 
Taking the square root of Imv an d Ih2R gives the Poisson errors in these counts, denoted 
ay and ctr, respectively 

The possible effects of using a different template spectrum other than the Arcturus 
atlas for the normalization procedure were explored using synthetic spectra. Model red 
giant spectra (normalized to their continua) were computed for effective temperatures of 
T e fj = 3800, 4130, 4540, and 5250 K, which span the range in temperature encompassed by 
the stars in our sample. Mismatches in spectral type associated with using these templates 
in the normalization of the observed spectra were found to introduce possible systematic 
effects in derived (V/R) values on the level of ~ 0.005, which is relatively small compared to 
uncertainties due to photon-statistics (see Table 1). For the random errors associated with 
the normalization process, which vary according to where the fitting points are set outside 
the H and K lines, the errors are typically 0.007 on (V/R)k and 0.003 on (V/R)h, which 
are again considerably smaller than Poisson noise in the photon counts. While the exact 
template used for normalization is not critical, it is important that the normalization be 
done. A test of the (V/R) values in the unnormalized spectra revealed differences of -0.02 
to -0.05 for the (V/R)h depending upon the color of the star. 

As noted, we have simultaneous spectra of the Ca II K emission line for each H-line 
observation listed in Table 1, these two features falling in different echelle orders of the same 
CCD exposure. The (V/R)k value of the Ca II K emission feature for each observation, and 
its error ex, are also reported in Table 1. These values were taken from Smith & Shetrone 
(2004) and Shetrone et al. (2008), with the exception of the June 2008 measurements. Several 
(V/R)k values were taken from Smith & Shetrone (2000) and adjusted with a zero point 
offset of +0.026 to correct for differences between normalization procedures, as described in 
Shetrone et al. (2008). The Smith & Shetrone (2000, 2004) papers did not include errors in 
the (V/R)k data. To estimate these we first determined an average difference between the 
errors in (V/R)h and (V/R)k from the Shetrone et al. (2008) data and the H analysis done 
here. This average difference, e# — ex, is 0.005. This offset was applied to the e# values 
derived from our analysis of the older Ca II H line spectra to estimate corresponding values 
of €k- These new e K error estimates, the errors from Shetrone et al. (2008), and the e# 
errors from this analysis are listed in Table 1. 

Photometry in the Johnson BV system was taken from the General Catalogue of Pho- 
tometric Data (GCPD; Mermilliod et al. 1997), accessed from an online catalog^. The B — V 
color of each star, as well as the absolute visual magnitude derived from Hipparcos parallaxes 
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(assuming zero interstellar absorption), are listed in Table 2. The parallaxes were tabulated 
in the earlier papers of Smith & Shetrone (2002, 2004). The distance range among most of 
the stars in our sample is 20-200 pc, with 9 stars of known parallax being beyond 200 pc. 
Also listed in Table 2 is the heliocentric radial velocity v r as obtained from the SIMBAD 
database and the metallicity from Mc William (1990) of each star. These properties are all 
consistent with our sample being Population I giants. 

In Figure 1 the asymmetry (V/R)h of the Ca II H line is plotted versus the B — V 
color of each star. In cases of stars observed multiple times, each observation is plotted 
separately, and consequently this diagram contains the effects of any time variability in the 
H line asymmetry. With the exception of one observation of HD 203387, only stars with 
B — V > 1.2 exhibit red-dominant asymmetry where (V/R) H < 1.0. The most extreme red 
asymmetries where (V/R) H < 0.7 are only seen in stars with B — V > 1.5. Figure 1 can 
be compared to Figure 2 from Shetrone et al. (2008), which gives an analogous plot for the 
asymmetry parameter of the Ca II K line. Both plots illustrate two groupings of data points 
spanning similar asymmetry and color ranges: one in the region 1.0 < (V/R) < 1.3 and 
0.85 <(B -V) < 1.05, and the other with 0.7 < (V/R) < 1.4 and 1.2 < (B - V) < 1.45. 
However, data for giants redder than B — V = 1.45 are unique to Figure 1, which represents 
a more diverse sample of stars than that in Shetrone et al. (2008). 



3.2. The Issue of Interstellar Absorption 

Stellar H and K emission line profiles can be affected by Ca II absorption in the in- 
terstellar medium (ISM). Although we do not have spectra of the interstellar Ca II lines in 
the directions of our program stars, several types of evidence indicate that the asymmetry 
trends seen in the (V/R) data are not attributable to systematic interstellar corruption of 
the stellar Ca II lines. 

If interstellar absorption was dominating the observed (V/R) values then the sense of 
the asymmetry should correlate with the velocity of each star relative to the ISM (as in 
Figure 5 of Bohm-Vitense 1981 for the Mg II /c2 emission). Several diagrams were prepared 
to test whether a similar correlation is present in the Ca II data. Plots were made of the 
average (V/R) K and average (V/R) H values of each star versus heliocentric radial velocity 
v r ; no trend is seen (Figure 2). Furthermore, a version of Figure 1 was made using different 
symbols to depict stars in different radial velocity ranges: v r < —20 km s _1 , —20 < v r < +20 
km s _1 , and v r > 20 km s -1 . There is no evidence that stars with positive or negative v r 
occupy systematically different regions of the diagram. We conclude that the Ca II emission 
lines are more resilient to ISM corruption than found by Bohm-Vitense (1981) for the Mg II 
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h and k emission. 

As a second test, the ISM absorption Ay towards most program stars was calculated 
using the numerical recipe of Hakkila et al. (1997), which takes into account both the 
distance to each star and the Galactocentric coordinates of each star. Many of the stars 
from Table 1 are within 75 pc of the Sun, and within a local bubble inside of which ISM 
extinction is considered to be small or negligible (e.g., Sfeir et al. 1999; Luck & Heiter 
2007). The equivalent width of the interstellar Ca II K line as seen against early- type stars 
within 100 pc of the Sun tends to be less than 10 mA (e.g., Vallerga et al. 1993), although 
exceptions occur. The values of Ay derived according to the Hakkila et al. (1997) model, 
modified such that Ay is set to zero within 75 pc, are listed in Table 2. For the majority 
of program stars, the estimated A v < 0.1 mag. Furthermore, those stars with the largest 
values of Ay do not occupy systematically discordant positions in the H and K line (V/R) 
versus (B — V) diagrams. In summary, considering both the radial velocities of the program 
stars, and the estimated interstellar extinction towards them, we see no systematics of the 
type that would be expected if ISM corruption of the Ca II H and K emission profiles was 
governing the observed (V/R) values. 



3.3. Differences Between the H and K Line Asymmetries 

Despite the similarities between Figure 1 of this paper and Figure 2 of Shetrone et 
al. (2008), there are nonetheless differences between (V/R)h and (V/R)k for many stars. 
Some illustrative examples of asymmetry differences are shown in Figures 3 and 4. The first 
of these shows spectra of HD 156283 obtained on 22 June 2008 using the 2.7-m telescope, 
wherein the (V/R)k ratio is 0.99 and (V/R)h is 1.07. At this time, although the K emission 
line was nearly symmetrical, the H line showed a blue-enhanced asymmetry in which the 
violet H 2 y peak was stronger than the red peak. Figure 4 shows a spectrum of HD 21552 
obtained on 3 October 2001 using the 2.1-m telescope. The Ca II (V/R) K ratio for this 
spectrum is 0.94 while (V/R) H = 1.01. In this case, although the H line is nearly symmetric, 
the K emission profile is stronger in the red peak. 

The difference between the Ca II H and K line V/R values was calculated for every 
observation of each star in the sample. The resulting set of difference values is plotted versus 
the H line (V/R) H in Figure 5; each observation of each star is plotted, so some stars are 
represented by multiple points. The error bars were calculated by summing en and ex in 
quadrature. The value of (V/R)h is larger than (V/R)k by an average offset of +0.039, 
which is shown as a dashed line in Figure 5. However, it can be seen that there are trends 
within this diagram, such that a single average offset is not a good fit to the entire data set. 
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The difference between the H and K line Vj R values are shown in Figure 6 versus B — V 
color. As with Figure 5, asymmetry values from individual spectra are plotted, so that some 
stars are depicted with multiple data points. The linear best fit to the data is plotted as 
a dashed line, which has a slope of 0.162 ± 0.017. However, this fit does not seem to be 
the best representation of the data. In order to smooth over possible intrinsic variability, 
an average of the (V/R)h — (V/R)k values was calculated for each star. These averages 
are shown plotted versus B — V in Figure 7, so that unlike Figure 6, each star is depicted 
by just one data point. In order to further smooth the data to look for trends, a weighted 
mean was made of every five points in Figure 7 binned according to B — V. A continuous 
line shows the locus of these binned averages. Rather than a linear relation, these binned 
data exhibit a quasi-step function with a reasonably abrupt change at (B — V) ~ 1.25. The 
two horizontal lines in the figure show the average value of the asymmetry difference among 
giants with colors on either side of this "divide." 

The color-magnitude diagram for our sample of red giants is plotted in Figure 8. The 
data points are coded according to the significance level of the (V/R)h — (V/R)k values; open 
circles and open squares represent stars for which the average asymmetries differences are 
((V/R) H - (V/R) K ) > la and ((V/R) H - (V/R) K ) < la, respectively, while filled triangles 
represent stars for which \ ({V/R)h — (V/R)k)\ < 1&- Red giants for which (V/R) is notably 
greater in the H line than the K line (open circles) tend to have colors of (B — V) > 1.2, 
while stars for which the reverse is the case (open boxes) are mainly of (B — V) < 1.2, 
although they do extend as red as (B — V) = 1.4. Giants with similar H and K line 
asymmetries (filled triangles) are evenly spread across the entire B — V range covered by 
the data sample. Thus, stars with similar asymmetry differences tend to group in specific 
regions of the color-magnitude diagram, albeit with some overlap, particularly in the color 
range 1.2 < (B — V) < 1.4. This is a further reflection of the trend seen in Figure 6. 



3.4. Photospheric Effects 

The Ca II H and K line chromospheric emission components are superimposed on a pho- 
tospheric absorption spectrum that varies with both wavelength and effective temperature. 
Could the differences between (V/R)u and (V/R)k be due to differences in the background 
photospheric spectrum in the vicinity of the H and K lines? We have investigated this ques- 
tion through the use of synthetic photospheric spectra calculated using the 2007 version of 
the program MOOG (Sneden 1973). These spectra were generated for a model atmosphere 
with T eff = 4250 K, logg = 1.75, and [Fe/H] = -0.1, a combination of parameters that is 
typical of the red giants in our sample. Synthetic photospheric spectra were generated for 
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three linelists: one with all of the photosphere lines, a second list from which the H and 
K lines had been removed, and a third that contained only the Ca II H and K lines. To 
simulate a symmetric double-peaked H or K emission feature, two equal-height Gaussian 
profiles were combined having separations in wavelength that are typical of the observed 
spectra. In accordance with the optically thin case, the relative strengths of these H and K 
emission profiles were set to be 1:2, with the peak flux in the K emission being 15% that 
of the continuum level of the photospheric spectrum on which it is superimposed. Mathe- 
matically, the flux of these artificial peaks was simply added to the synthetic photospheric 
spectrum computed with the full line list. 

A decomposition of the resultant model spectrum is presented in the various panels 
of Figure 9, which show from top to bottom: the H and K photospheric lines alone, all 
photospheric lines other than H and K, the complete photospheric spectrum, the artificial 
H and K emission profiles, and the full photospheric spectrum plus the artificial emission 
lines. There are numerous weak and moderate strength absorption lines throughout the 
wavelength range 3925 < A < 3980 A, but once they get superimposed on the very broad 
and deep photospheric profiles of the H and K lines, the amplitudes of these weaker lines 
are greatly reduced. Since the total photospheric flux in the cores of these lines is relatively 
small, the effect of the photospheric background on the chromospheric emission profiles is 
very modest. In Figure 9, the asymmetry induced by this effect is barely discernible. 

To complement the simulation, photospheric models were also computed for a range 
of effective temperatures. As with Figure 9, simulated H and K emission lines, both with 
iy/R) = 1 and a K line twice as strong as the H line, were added to the set of synthetic 
photosphere spectra. From these superimposed spectra the values of (V/R)k, (V/R)h, and 
(V/R)h — {V/R)k were calculated. In addition to simulations with symmetric emission 
profiles, calculations were also done for input values of (V/R) h ,k = 1-2 and 0.6, as well 
as a range of wavelength separations between the two maxima in each bimodal profile. In 
all cases the input H and K emission profiles had identical asymmetry, such that initially 
(V/R)h — {V/R)k = 0. The value of this difference after superimposing the model photo- 
sphere spectrum illustrates the extent to which photospheric corruption might be affecting 
the observed asymmetry differences. The resultant values of {V/R)h — {V/R)k are plotted 
versus B — V in Figure 10. Modest systematic variations are seen as a function of (B — V), 
which are largely driven by photospheric corruption of the K line rather than the H line. 
However, the magnitude of these asymmetry differences is much smaller (by a factor of 5) 
than the star-to-star differences, and their trend with (B — V), than we actually observe. 
Although this modeling approach is no substitute for a self-consistent treatment of radiative 
transfer through a combined model photosphere plus chromosphere, it does suggest that the 
iy/R) asymmetry variations among red giants are not an artifact of photospheric differences, 



- 9- 



but reflect physical differences in the stellar chromospheres. 

In computing the synthetic spectra it was assumed that there is no systematic velocity 
shift between the simulated emission profiles and the photospheric spectrum. As noted 
above, in the data reduction procedure all of the observed spectra were shifted so that the 
photospheric lines are at the rest wavelengths as defined by the Arcturus atlas. The average 
central wavelengths of the H and K emission profiles were measured in this rest frame for 
the 8 stars observed at the 2.7-m telescope in 2008. There was no systematic displacement 
of the emission line centers from the photospheric reference frame. The average deviation 
was found to be —0.02 km s" 1 with a standard deviation of 1.0 km s _1 (1 pixel = 2.4 km 
s _1 ), and there is no trend with (B — V) color. 

4. Conclusions 

Values of the asymmetry parameter (V/R)h for the Ca II H chromospheric emission 
line have been measured for a large sample of red giant stars. The results are presented in 
Table 1, plotted as a function of B — V in Figure 1, and in a series of other diagrams are 
compared to V/ R values of the Ca II K line based on the same spectroscopic data set (see 
Shetrone et al. 2008). Many of the spectra exhibit opposite asymmetries between the H and 
K line emission peaks. The H line (V/R) H value is larger than the (V/R) K ratio for the 
majority of stars in our sample, however, this difference shows a trend with B — V color, 
being largest among giants with (B — V) > 1.25. 

The observations can be interpreted within a framework in which the H and K emission 
lines are formed at different levels in a stellar chromosphere, thereby allowing for a compar- 
ison between the motions of these layers. Asymmetry in either the H or K line, i.e., V/R 
different from unity, is generally attributed to Doppler shifts resulting from a differential 
velocity field in the chromosphere. In this Doppler picture, gas that produces the H 2 y and 
H 2 r peaks is taken to be moving with a different velocity to the gas that contributes to the 
central H 3 reversal. This causes the H 3 reversal to be redshifted or blueshifted with respect 
to the H2 peaks, resulting in an emission line asymmetry. Such a scenario has been invoked, 
for example, to interpret H and K line asymmetry differences among dwarf stars (Rauscher 
& Marcy 2006). 

On the Sun the V/R values for the H emission line follow a slightly skewed distribution 
with a maximum at unity, a minimum of (V/R)h ~ 0.5, and an extended tail to values 
greater than 1.7 (Cram & Dame 1983). There is only a weak correlation between (V/R)h 
and the strength of the H emission, from which Cram & Dame infer that "the V/R ratio 
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of the integrated solar H line is not due specifically to bright (e.g., network) elements". 
Furthermore, the H 3 reversal is found to be consistently displaced away from whichever of 
the H 2 y and H 2 ^ peaks is brightest, which indicates that Doppler velocity shifts of the H 3 
reversal are an important factor governing the asymmetry of the H emission profile (Cram 
& Dame 1983). Enhanced emission in the H 2 y peak appears to be associated with shock 
heating (e.g., Beck et al. 2008). In regions of the quiet Sun, the H 2 y and H 2 # emission have 
even been observed to disappear, indicative of cool elements within cells of the chromospheric 
network (Rezaei et al. 2008). 

Among the majority of dwarf stars with spectral types in the range K7-M5, Rauscher & 
Marcy (2006) show that V/R > 1 in both the H and K lines, with a preference for (V/R)k > 
(V/R) H . In addition, the K 3 reversal tends to have a systematically greater redshift than the 
H 3 feature. It is worth quoting the conclusions of Rauscher & Marcy on this point: "While 
the unequal degeneracies in the upper levels producing these lines cause somewhat differing 
optical depths, these data emphasize that the H and K lines are produced in slightly different 
regions of the chromosphere and are therefore governed by slightly different line transfer 
and kinematics." The redshift of the central reversals is attributed by them to micro-flare 
activity that heats chromospheric gas and propels it outwards. As it moves upward it cools 
and decelerates, resulting in the gas that contributes to the central reversal being redshifted 
relative to the faster moving, hotter gas below it that produces the emission peaks. The 
velocity differential between these regions of gas is on the order of several tenths of a km 
s -1 . The greater redshift of the K 3 reversal compared to the H 3 reversal is consistent with 
the former arising from slightly higher altitudes in the chromosphere, where deceleration has 
proceeded to a greater degree. 

Clues to the systematics of the chromospheric motions at work within our sample of 
red giants can be found in Figures 11 and 12. Average values of (V/R) H and (V/R) K were 
calculated for each star using the data in Table 1. The main panel of Figure 11 shows these 
average asymmetries plotted versus one another. In Figure 12 both ((V/R)h) and ((V/R)k) 
are plotted versus the average value of the asymmetry difference ((V/R)h — (V/R)k)- This 
diagram is analogous to Figure 5, except that all data for a given star are averaged, so that 
each giant is represented by a single point. These figures map out the relative behavior of 
the H and K emission line profiles. 

Among giants with the largest values of (V/R)h (> 1-1) this asymmetry parameter is 
nearly equal in value to (V/R)k', such stars exhibit a blue-dominant asymmetry in both the 
H and K emission profiles, and scatter closely around the locus {V/R)h = (V/R)k which is 
shown by the straight line in Figure 11. Among such stars it seems that the H and K lines 
are evincing similar velocity fields and that there is no evidence of systematic outflow in the 
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chromosphere. 

Turning to giants with 1.0 < (V/R)h < 1-1, there is a noticeably larger range of 
K line asymmetries and one starts to see a systematic shift of the K line profiles to red- 
dominant asymmetries, i.e., to (V/R)k < 1-0. A significant fraction of giants with 1.0 < 
(V/R)h < 1-1 also have (V/R)k < 1-0, and so exhibit opposite H and K emission profile 
asymmetries. Among such stars we infer that an accelerating outflow has become established 
within a region of the chromosphere that spans altitudes at which the K 2 and K 3 features 
are formed. We interpret the later feature as being formed at higher levels, in regions of 
faster outflow, than K 2 . This causes the central absorption reversal to be blueshifted with 
respect to the emission profile, thereby giving a characteristic asymmetry of (V/R)k < 1-0. 
However, the outflow does not yet extend to the deeper altitudes where the H 3 feature is 
formed. It is therefore a change in the K line asymmetry that drives the non-zero values of 
((V/R)h — (V/R)k) among giants having 1.0 < (V / ' R)h < 1-1. As a result, many of these 
stars fall below the straight line in Figure 11. This phenomenon is reflected by a plateau of 
approximately constant (V/R)h ~ 1.1 ± 0.05 in the upper panel of Figure 12, whereas the 
same stars fall along an inclined locus in the lower panel. 

Eventually the accelerating outflow extends downwards in altitude to regions where the 
H3 feature is formed, resulting in (V/R)h now decreasing to a value below unity. This 
brings about a maximum in the mean difference ((V/R)h — (V/R)k) of ~ 0.2. Both (V/R)h 
and (V/R)k are now less than unity and we infer that an accelerating wind now extends 
throughout the chromospheric regions where the K 3 , K 2 , H 3 , and possibly the H 2 features 
originate. As the value of (V/R)h drops below 0.8 the difference in V/R between the H and 
K lines diminishes, and for stars with (V/R)h ~ 0.4-0.5 the asymmetry parameters for the 
two lines become approximately equal. Such giants may be those with the largest velocity 
gradients in the chromosphere. They have returned to the vicinity of the straight line plotted 
in Figure 11. 

Figure 1 shows that the systematic trends outlined in the previous paragraphs corre- 
spond approximately to a sequence in B — V color. Thus it may be that as stars evolve 
up the red giant branch they travel along an evolutionary sequence in Figures 11 and 12, 
with winds first originating at high altitudes in the chromosphere, and then becoming more 
established and extending to lower altitudes. 

However, it is also known that V/R for the H and K lines can be time variable (Shetrone 
et al. 2008), and this is likely to impose scatter on any long-term evolutionary sequence. An 
excellent example is the star Arcturus, which both Gray (1980) and Qin & Li (1987) find to 
be variable over timescales of hours to months. This indicates that red giant atmospheres are 
likely subject to both short and long term dynamical modification. Changes on a timescale 
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of hours were interpreted by Qin & Li (1987) as being due to a "chromospheric eruption" 
that may have ejected mass rapidly through the chromosphere and possibly to escape from 
the star. 

In the solar chromosphere significant differences in the ratio of H and K line source 
functions can occur over a range in height of approximately 500-700 km (Linsky 1970). 
Comparing to a G giant of solar effective temperature and mass, and radius R = 40R Q , the 
chromospheric scale height might be expected to vary roughly as the inverse surface gravity, 
i.e., as R 2 . Thus, variations in the ratio of the H and K line source functions might extend 
over an altitude range of ~ 800, 000 km. A disturbance moving at 5 km s _1 could traverse 
such a region in about 44 hours. Although we do have multiple observations of many stars, 
the typical time difference between consecutive spectra is months or years rather than days. 
Thus, in general, we don't have fine enough temporal coverage to detect possible systematic 
variations in H or K line asymmetries due to passage of a short-term "event" through the 
chromosphere. In the case of HD 27697, the two spectra obtained in 1998 on consecutive 
nights reveal no evidence for any significant differences in (V/R) for either the H2 or K 2 
features. Similarly, spectra of HD 82635 taken 18 nights apart in 2007 show no significant 
evidence of variability. In 2006, the data show that both the H 2 and K 2 {V/R) values for 
HD 222404 differed significantly between 1 September and 4 October, but our observations 
were a month, rather than several days, apart. 

It should be noted, however, that the H and K line observations of Arcturus by Qin & Li 
(1987) may depict a very energetic version of a short-timescale phenomenon. The asymmetry 
indices measured by them (using their notation where r refers to the intensity in an emission 
line peak) vary within the range 0.87 < r(H 2 y)/V(H 2r ) < 1.21 and 0.79 < r(K 2 ^)/r(K 2r ) < 
0.93 over a period of 4 hours, and all their observations show r(H 2 y)/r(H 2r ) > r(K 2 y)/r(K 2r ). 
Furthermore, in 4 of their 6 spectra r(K 2 y)/r(K 2r ) < 1.0 while r(H 2 y)/r(H 2r ) > 1.0, i.e., the 
H and K emission lines evince opposite asymmetries. On the basis of this example, it follows 
that short-term variability could indeed contribute to the star-to-star scatter in asymmetries 
seen in the McDonald spectra. Such effects may be superimposed on a much longer-term 
evolutionary sequence that produces the trends with (B — V) color seen in Figures 11 and 
12, and which may be associated with the onset of pulsation and/or the driving of a large 
scale wind. 

This project was completed as part of the McDonald Observatory REU program and was 
supported under NSF AST-0649128. This research has made use of the Simbad database, 
operated at CDS, Strasbourg, France. We would like to thank E. Luck for providing the 
extinction code, hakkila, used in this work. 
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a Metallicities taken from McWilliam (1990). 



-25 - 



1.4 



1.2 



g 1 
> 

QJ 0.8 



0.6 



0.4 




0.6 



0.8 



1 1.2 

B-V color 



1.4 



1.6 



Fig. 1. — The value of the Ca II H line from Table 1 is plotted versus the B — V 
color (Table 2). The single star symbol represents Arcturus. Where multiple observations 
are available for a star, the result for each individual spectrum is shown as a separate data 
point. Thus a number of stars are represented by multiple points. 
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Fig. 2. — The H and K line asymmetry parameters (V/R)k and (V/R)h versus heliocentric 
radial velocity. The lack of correlations in this figure is presented as evidence that the 
asymmetry trends seen in this, and previous papers, are not dominated by the effects of 
interstellar absorption. 
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Fig. 3. — The order of the normalized spectrum of HD 156283 that contains the Ca II H and 
K lines is plotted in the upper panel. This spectrum was observed on 22 June 2008 using the 
2.7 m telescope at McDonald Observatory. In the bottom left window, an expanded version 
of the core of the K line is plotted, to clearly show the near symmetry of the emission peaks. 
The core of the H line is plotted in the bottom right window, and reveals the blueward 
emission peak to be stronger than the redward peak. 
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Fig. 4. — The echelle order of the normalized spectrum of HD 21552 containing the Ca II 
K line is plotted in the top panel, with the H-line order shown in the bottom window. 
This spectrum was observed on 3 October 2001 using the 2.1 m telescope at McDonald 
Observatory. Whereas the H emission line is nearly symmetric, the K line shows a slightly 
red-enhanced emission profile. 
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Fig. 5. — The difference between the Ca II H line and K line V/R values plotted versus the 
measured H line V/R. The data are from Table 1, so certain stars may be represented by 
multiple points in this diagram. The single star symbol represents Arcturus. The dashed 
line shows the mean value of the V/R parameter difference. 
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Fig. 6. — Observed values of (V/R) H — (V/R) K versus B — V color. The data are from Table 
1, so some stars show as multiple points in this diagram. The single star symbol represents 
Arcturus. The dashed line shows a least squares linear fit. 
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Fig. 7. — The average (V/R)h — (V/R)k value for each star listed in Table 2 versus B — V 
color. A continuous line connects the weighted mean of every 5 points binned according to 
B — V. This curve shows a rather abrupt change at B — V ~ 1.25. The mean values on 
either side of this "break" are shown by the horizontal lines. 
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Fig. 8. — The color magnitude diagram of the present sample of giants listed in Table 2. 
The open circles and open squares represent data points for which the average {(V/R)h — 
(V/R)k) > Iff or < —la - respectively. Filled triangles represent data points for which 
\((V/R) H -(V/R) K )\<la. 
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Fig. 9. — Synthetic spectra at wavelengths encompassing the Ca II H and K lines for a star 
with T e fj = 4250 K, logg = 1.75, and [Fe/H] = -0.1. Horizontally paired panels correspond to 
the regions of the K (left) and H (right) lines. Running from top to bottom the panels show: 
(i) the H and K photospheric line profiles alone, (ii) a photospheric spectrum based on a line 
list that excludes H and K, (iii) the full photospheric spectrum, (iv) simulated chromospheric 
emission lines with a ratio of 2:1 in the K:H line strengths, and (v) the summation of the 
complete synthetic photospheric spectrum plus the simulated emission components. 
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Fig. 10. — The asymmetry difference (V/R)h — (V/R)k, found after taking H and K emission 
lines of equal initial (V/R) and superimposing them on a photospheric spectrum, versus the 
(B — V) color of the photosphere. Symbols denote the morphology of the intrinsic H and K 
emission lines, which are parameterized by their initial (V/R) values and the peak-to-peak 
separation of the bimodal profiles. 
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Fig. 11. — The large upper panel shows the average value of (V/R)k versus the average 
(V/R)h for each star listed in Table 2. The lower panel plots the average value of the 
asymmetry difference (V/R)h — (V/R)k versus average (V/R)h- The solid lines are the 
locus for (V/R) H = (V/R) K . 



-36- 




Fig. 12. — The average value of (V/R)h and the average (V/R)k for each star in Table 
2 plotted versus the average value of the asymmetry difference (V/R)h — (V/R)k- Unlike 
Figure 5, each star is represented by a single point in this figure. 



